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FOREWORD 


The investigation conducted at USC concerning the roie of the large scale, 
quas i -ordered structures of a turbulent jet In the generation of noise has been 
jointly supported during the period of this contract by the Department of 
Transportation under grant number DOT-OS-00002 and the National Aeronautics and 
Space Administration. The research conducted under NASA Contract NAS3“1?857 
and reported herein emphasizes those aspects of the problem that are related to 
the near rather than far field. The structure and content of this report re- 
flects this fact. 
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SUMMARY 


Previous flow visualizations of turbulent jets at low Reynolds numbers 
have clearly indicated the presence of vortex ring-like structures and the 
nature of their interaction. 

A series of measurements of near field pressures and turbulent velocity 
fluctuations were made in a Jet having a Reynolds number of about 50,000 in 
order to investigate more quantitatively the character and behaviour of these 
structures and to ascertain their importance to the Jet noise problem. 

Auto- and cross-correlation measurements of the pressure fluctuations 
Just outside of the Jet yielded information about the average separation dis- 
tance between neighboring structures and their convection velocity. 

Coalescence between neighboring structures near the nozzle was detected 
by auto-correlat Ion and two-point eduction techniques, however, its detection 
further downstream proved to be unsuccessful. As a consequence, a completely 
new statistical technique has been initiated in order to study this problem. 

It was found that the process of interaction between vortices can be in- 
hibited by artificially exciting the shear layers with periodic disturbances, 
of certain frequency. The turbulent fluctuation amplitudes measured at four 
diameters downstream decreased considerably. 

Finally, it was observed that the passage frequency of the structures 
decreased with x in a similar manner as the frequency corresponding to the 
maximum intensity radiation emanating from the same value of x. 

The work is incomplete in the sense that whereas some new information was 
found about the large scale structures of jet turbulence, little quantitative 
information about their spatial character and their mutual Interaction could 
be obtained; furthermore, as yet no convincing direct connection between far 
field sound and the large scale structures in the flow could be demonstrated. 


INTRODUCTION 


For the past two decades, Jet noise research has been dominated by the 
work of lighthMl (ref. 1). In his theory the Jet is modeled as a passive 
fluid medium with sound speed a^ , containing a distribution of "equivalent" 
acoustic sources whose strength is equal to the scalar contraction of the 
local fluid stresses present in the Jet. The density field, , must then 
satisfy the acoustic wave equation (ref. 1): 




- 
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In Llghthill's model the fluid stress, Tji , is considered to be distributed 
over a large number of random, uncorrelatea “eddies". 


There is a growing body of evidence indicating the existence of large 
organized structures within turbulent shear flows. Most recently in a Ph.D. 
thesis at the Rhein-Westfal ia University, Aachen, West Germany, Gunther Neuwerth 
has demonstrated the presence of large coherent structures in a jet at high 
subsonic Mach numbers using an ingenious multispark Schl ieren system (private 
communication). Although it is still an open question whether or not these 
structures contribute directly to the far field noise, one can make reasonable 
arguments that they might. 

In an important recent work Crow and Champagne (ref. 2) detected large 
scale organized structures in a jet up to a Reynolds number of about 80,000 
by visualization techniques. Furthermore, they used acoustically driven stag- 
nation chamber resonances to excite the jet (Re = 100,000). They measured the 
jet response by placing a hot wire anemometer within the irrotational core where 
velocity perturbations would, presumably, reflect gross features of the local 
shear layer turbulence. They observed that the natural peak frequency occurred 
at a Strouhal number (based on jet exit velocity and exit diameter) of 0.3 and 
near the end of the potential core (x/D = 4). However, when the Jet was excited 
at the first harmonic (St = .6) , the fundamental was amplified at a still 
greater rate. 


On the basis of their results. Crow and Champagne tended to regard the jet 
flow in terms of a column instability. In fact, in an unpublished work Crow 
formulated a “wave antenna model" in place of Lighthill's compact, random source 
model, in order to estimate the far field noise of a turbulent jet. 


However, the presence of the large scale structures might be interpreted in 
terms other than that of wave instabilities. The measurements of Lau, Fisher 
and Fuchs (ref. 3) indicate that the shear layer is made up of a “row" of dis- 
crete rings of vorticity. This is consistent with the visualization of Winant 
and Browand (ref. 4) made in a two-dimensional free mixing layer. They found 
that due to the linear wave instability present near the origin of the shear 
layer, the initial vortex shear layer rolls up into cylindrical lumps of concen- 
trated vorticity and the adjacent lump? roll around each other and coalesce into 
a single large lump. This “pairing" process was postulated as the major mechanism 
of turbulent shear layer growth. Vortex pairing has also been observed by Rock- 
well and Niccolls (ref. 5) in a two-dimensional jet using low Reynolds number 
flow visualization techniques and by Browand in an axisymmetric water jet as 
reported by Laufer (ref. 6). 


There are well justifiable reasons for investigating vortex ring coalescence 
as a possible mechanism of far field noise production. Powell has shown (ref. 7) 
that when a free vortex ring stretches, it radiates acoustically like a dipole 
and the far field particle velocities are proportional to the second time deri- 
vative of its Kelvin impulse. During a pairing between two rings, the absence 
of an external force demands the Kelvin impulse of the pair is invariant. 

However, their far field, because of the Stokes effect, degenerates into a 
quad rupole- 1 i ke behaviour. 
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On the basis of the preceeding discussion it is apparent that the question 
of the large scale structures as potential "acoustic" sources is a relevant one 
indeed. The principal aim of the research described in the present report is 
directed toward a detailed study of this question. For this purpose it is 
necessary to develop techniques for the detection of the organized structures 
inside the jet while simultaneously measuring their contributions to the far 
field noise. One should remark at this point that conventional cross correlation 
techniques involving one transducer in the flow and one in the far field are not 
particularly helpful, since single point flow measurements cannot provide suf- 
ficient information about the large structures. 

The experimental problem has been approached in three parts. First, the 
behaviour of the large structures has been studied with particular emphasis on 
how their presence is related to the near field pressure and velocity signatures. 
In connection with this study the effect of artifically introduced disturbances 
near the initial jet development has also been included. Second, specific 
features of the radiated noise using a reflector microphone technique have been 
explored (primarily under DOT sponsorship). Finally, a search for a technique 
to relate the large scale structures to the noise field was undertaken. 

In the first part of the Results and Discussion section, it is shown that 
the flow visualization techniques work well in detecting the behaviour of the 
large scale coherent structures In low Reynolds number jets. However, diffi- 
culties arise for moderate and high Reynolds numbers and the following two sub- 
sections describe the various approaches and techniques applied in an attempt 
to bring to light the characteristic features of the coherent structures by 
means of near field measurements. The subsequent discussion includes results 
obtained by the application of initial disturbances artificially introduced 
near the nozzle exit. Finally, measurements obtained primarily with a reflector 
microphone technique are presented with particular emphasis on how these meas- 
urements might relate to the large scale structures. 


FACILITY AND INSTRUMENTATION 


Model Jet 


Most of the experiments described in this report have been carried out 
in a low speed jet facility. It consists of a vacuum cleaner blower, a settling 
chamber with damping screens and a contraction section with a circular nozzle. 
The contraction ratio is 36:1 and the nozzle diameter is 2.5^ cm. The 
attainable maximum jet velocity is 70 m/sec, but most of the work was carried 
out at velocities near 30 m/sec. This corresponds to a Reynolds number based 
on exit velocity and nozzle diameter of 5.2 x 10^. 

The exit boundary layer is estimated to be less than .05 diameter. The 
exit mean velocity was found to be flat to better than 0.5% and the exit plane 
turbulence level was measured at 0.3%. This relatively high level was consid- 
ered satisfactory for the experiments planned In the facility. 
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Anechoic Chamber 


Far field measurements quoted in this report were carried out in the USC 
Anechoic Chamber facility under the sponsorship of the Department of Transpor- 
tation (ref. 8). Briefly, the facility is of a blowdown type consisting of 1 
or 2 inch jet nozzles with carefully designed valving and settling chambers 
exiting into a large anechoic chamber. Additionally, provisions have been 
included for extensive monitoring and measurement of the parameters which have 
been of importance for characterizing the properties of the turbulent jet. 


Directional Microphone System 


The directional microphone system is based on the imaging property of a 
spherical reflector. The reflector used has a radius of curvature of 53" and 
an aperature of 36". The dimensions were chosen after a parametric study of 
the imaging properties of spherical reflectors. A 1/8" B S- K condenser micro- 
phone is used to pick up the focused acoustic waves at the geometrical image 
point. 

The reflector disc is made of aluminum plate backed by another aluminum 
disc with rubber like material sandwiched in between to damp out the natural 
vibration of the primary disc. A gun-sight telescope is mounted at the back 
through the center of the reflector for alignment purpose. 

Fig. 1 shows a schematic diagram of the traversing mechanism and the dir- 
ectional microphone system. This mechanism enables the directional microphone 
to scan automatically along the jet axis from a fixed location in the far field. 
At prescribed directions, the mechanism will pause for a few seconds to allow 
data to be recorded on tape before changing direction. The system also provides 
automatic refocusing at each prescribed direction. The scanning and refocusing 
are governed by a special electronic command system designed and built in our 
department . 


I nstrumentation 


Mean velocity measurements were made with a pitot tube and traversing mech- 
anism. Pitot tube total pressure was recorded off a capacitance-type differential 
pressure transducer. 

Near field pressure measurements were made using 1/8" capacitor microphones 
with a nominal sensitivity of -80 d8 re. iV'/ ^bar . Microphone output was 
preamplified by a cathode follower. 

Shear layer measurements were made with constant temperature hot wire 
anemometers. Typically the overheat ratio was 2, the wire was ,0001 inch dia- 
meter Platinum, and the aspect-ratio was 400. In the presence of a large ratio 





of mean velocity to r.m.s. fluctuation velocity, the wire responds to the velo- 
city component parallel to the mean flow. Consequently, a single wire was 
Sufficient to record axial velocity fluctuations. A cross wire arrangement was 
used when it was necessary to decompose the turbulence Intq axial and radial 
components. In all cases King's law, == A + B ( u cos 0)^ was assumed for cali- 
bration purposes, where E (t) is the hot wire signal; u(t) Is the flow speed; 0(t) is the 
angle between the velocity and the normal to hot wire and A, B are constants of cal 
ibration. In the case of the crossed wire probe it is necessary to add and subtract 
the outputs of the two wires in a prescribed manner to derive the orthogonal 
velocity components. 

Most of the statistical data processing was performed digitally using both 
hard wired equipment and a general purpose computer. Correlations and simple 
signal eductions were accomplished on a digital correlator, an instrument which 
had the capability of performing both auto-and cross-correlations, as well as 
signal recovery used for conditionally sampled time averages. An oscilloscope 
was adjusted to trigger on the slope and level of a reference signal, and the 
"gate out" was used to drive the "external trigger" on the correlator; for 
each trigger the correlator would average one frame of the data. In this way 
superimposed randomness would average to zero in favor of any wave form with a 
constant phase relative to the trigger condition. More sophisticated statis- 
tical analysis, such as joint probability distributions, was accomplished using 
the general purpose digital computer. Raw signal and calibration data were 
converted from continuous analog records to periodically sampled digitized 
records using an analog to digital converter subsystem of our own design. 


MEASUREMENTS AND TECHNIQUES 
Velocity 


Mean velocities were measured with the Pitot tube apparatus using standard 
techniques while fluctuation measurements were carried out with the constant 
temperature hot wires discussed in the previous section. 


Near Field Pressure 


There is always some question as to what a microphone measures in the 
presence of mean flow. Work by Siddon (ref. 9) and Fuchs (ref. 10) have indi- 
cated that there are flow configurations where microphones can record the static 
pressure fluctuations with sufficient accuracy. The microphones were placed 
just outside the edge of the jet where mean velocities were small, and were 
oriented so that the plane of the diaphragm was roughly parallel to the mean 
flow. 


Auto-correlations and cross-correlations of pressure signals were 
measured using the digital correlator. The microphones were located on 
the surface of a 10*^ cone intersecting the inside diameter of the 
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nozzle exit. Auto-correlations are shown at various downstream stations in 
Fig. 2. The r.m.s. level at each station is normalized to unity. Cross-cor- 
relations were measured with two microphones displaced a known distance, A 
displacement of 0.5 cm was selected for axial cross-correlations (Fig. 3) 
while azimuthal cross-correlations were measured along a cone angle of 10 
(Fig. 4). 


Special Statistical Techniques 

Experience indicates that the unambiguous detection of large, coherent 
structures Inside the shear layer is a difficult problem in signal recovery, 
especially downstream of the potential core. Time averaged techniques have 
proven successful only in detecting the first pairing; beyond that we have 
applied various eduction schemes developed recently at USC in connection with 
boundary layer and free shear layer research. 

Two Point Eduction .- In order to measure phase relationships internal to 
a coherent structure at various depths in the shear layer, two point eduction 
was used. This technique, identical to that used by Winant and Browand (ref. 
4), consists of sampling the properties of a large number of events occurring 
at a fixed position and at known time intervals before or after the occurrence 
of some judiciously chosen trigger conditions, measured at a different station 
and characteristic of the organized structure. Trigger position was inside 
the potential core, but close to the shear layer. The trigger condition 
selected was a positive signal slope and a signal level corresponding to one 
r.m.s. of the local fluctuation above the mean. This condition is a standard 
first attempt in conditional sampling of turbulent shear flows (ref. 23). At 
each station educted waveforms were recovered at the same downstream position 
but at varying radial positions. Educted waveforms appear in Fig. 5. Each 
figure corresponds to a particular downstream station, and in each case the 
indicated radial positions are relative, with positive positions in the dir- 
ection of the entrainment region. 

Joint Probability Measurements .- Although it is not possible to discern 
large scale turbulent patterns from single point measurements, it is possible 
to detect deterministic "subprocesses" in the turbulence from the Joint sta- 
tistics of orthogonal velocity components. The measurements were taken at the 
station x/D = 2, where the shear layer was fully turbulent, but where the 
vortex ring core was still reasonably localized so that there was room for 
vortex spreading. _Two radial positions were chosen: r/D = .49 and .56 

corresponding to u(r)/Ug = .79 and .60, respectively. 

The crossed wire probe was used to decompose the turbulence into radial 
and axial components. The hot wire signals were recorded on magnetic tape and 
then converted to a digital record for processing. From recorded calibration 
data, the digital records were linearized and stored as parallel records of 
radial (uj.) and axial (u^) velocities. Velocity space was divided into 
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2500 cells and the fraction of the total time the record occupied a given cell 
(u U|.) is the joint probability distribution P(Uj^, u^) evaluated at the 
cell and normalized by the peak value of the distribution. The resulting Joint 
probability distributions are shown In Fig. 6. The closed curves are contours 
o_f eq]ja1 probability. The origin of the coordinate axes is the mean velocity 
(uj^, U|.) . Generally speaking, the distributions consist of a single peak with 
a ridge, indicated as a dashed line, connecting the first, arid third quadrants. 


Time averages of the turbulence can be computed as moments of the Joint 
probability distribution. Specifically mean velocities are defined 



variances are defined 




Note that the probability distribution is normalized so that 

ff - / 

— oc? 

and hence, the correlation coefficient is defined 


( 1 ) 


( 2 ) 


( 3 ) 


( 4 ) 








'a' 


( 5 ) 


D i sturbance Generator 


The experiments of Crow and Champagne (ref, 2) show that gross properties 
of the Jet, such as axial velocity distribution, can be influenced by control- 
ling the Initial conditions of the jet. Their forcing technique, (acoustically 
driven stagnation chamber resonance) changes the shear layer structure directly 
since periodic exit plane modulation fluctuations induce a periodic series of 
ring vortices, a technique typically requiring forcing levels in excess of 

0.01 Ug. 

An early attempt was made to thermally drive the model air jet shear layer 
with an axisymmetric glow discharge. Thrs approach was subsequently abandoned 
due to the great difficulty in creating a uniform flow at atmospheric pressures. 
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More success was attained in driving the jet acoustically. A hemispherical 
chamber was fitted outside the nozzle contour to induce disturbances in the 
shear layer by acoustic pressure perturbations. Two speakers were used to 
excite a pressure field in the chamber, which communicated to the shear layer 
via a narrow gap, as shown in Fig. 7. In this way, an attempt was made to 
control the shear layer structures indirectly by providing an annular shaped 
acoustic wavefront that surrounds the entrainment region of the jet and where 
wave number j< is parallel to the jet axis. It was hoped that by driving 
the jet in this way, at the right frequency, a given vortex ring pairing could 
be phase locked with the acoustics, facilitating detection. This method of 
forcing is expected to differ' quai itatively from driving the jet core directly. 


RESULTS AND DISCUSSION 


Nature of the Turbu 1 ent Jet 


The model Jet exhibits general properties that are In good agreement with 
previous results. The centerline mean velocity decay is shown in Fig. 8 and 
is compared to the recent data of Crow and Champagne (ref. 2). Furthermore the 
linear spreading rate of the jet in the initial stage of development is depicted 
in Fig. 9» while Fig. 10 gives two typical mean velocity profiles in this region. 

The non-steady aspect of the flow field is, of course, the primary Interest 
here. The conventional root-mean-square velocity fluctuations will be described 
in a subsequent section; the presence of randomly occurring vortical structures 
will be discussed first. The most dramatic evidence of the existence of these 
vortical structures comes from flow visualizations. 

Two new, unpublished sequences of flow visualizations have been made in the 
laboratories at USC by Dr. F. K. Browand. The first of these is a film, in which 
the initial vorticity of the Jet is marked with dye (in a manner similar to that 
used in ref. 4) and its agglomeration into turbulent vortex rings is observable. 
The pairing process of Ref. 4 is easily observable In this film In the form of 
interactions of neighboring rings and can be seen as Fig. 2 of ref. 6. The 
second series of observations involved hydrogen bubble tracers and show sub- 
stantially the same phenomena, although with the added advantage of being capa- 
ble of providing more quantitative information. 

These visualizations of the unsteady turbulent structure of the jet motivated 
further research at higher Reynolds numbers. 


The Coherent Structures and Their Near Field 


Reynolds Number Effect .- As discussed in the previous section, visual techniques 
provide a convincing picture of the existence of the large scale, vortical 
structures and the process of their coalescence. At higher Reynolds numbers 
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(Re ISjOOP) visual techniques are much more difficult to use, and do not 
lend themselves easily for quantitative study. One might rely on the prin- 
ciple of Reynolds similarity in arguing that the general structure of the jet 
is not altered (the jet spreading angle is independent of the Reynolds number); 
the argument does not assure one that the shape and behavior of the large 
structures is the same for higher Reynolds number values. For this reason it 
is essential that one should develop alternate techniques for studying the large 
scale motion. It would be desirable to carry out these studies under high speed 
flow conditions where noise generation is readily detectable, keeping in mind 
that one's primary interest is establishing a connection between the large 
structures and the radiated noise. Because of the well known difficulties with 
using hot wire techniques at high speed, a compromise was made and a Reynolds 
number of about 50,000 was chosen for most of the reported measurements. It was 
noted namely that over a large Reynolds number range and at least down to values 
of 90,000 the normalized power spectral density of the overall radiated noise 
exhibits a similar behavior (Fig. 11). Although there is an observable shift 
of acoustic power from lower to higher frequencies with increasing Reynolds 
number, the location of the peak power is remarkably insensitive to Reynolds 
numbers under widely different experimental conditions. With the conjecture 
that the peak frequency is related to the time scale appropriate to the behaviour 
of an organized turbulent structure, one would expect to find the same type of 
structures In high and low Reynolds number jets. Admittedly this argument can- 
not be considered a conclusive one; nevertheless, if measurements at Re = 

50,000 show a picture of the large scale structure consistent with that 
obtained at Re = 5.000, this fact gives one more confidence to accept the 
Reynolds similarity argument for still higher values of the Reynolds number. 

Nature of the Large Structures .- It is a well accepted fact, based on 
past experience, that the understanding of conventional statistical measurements 
in a turbulent flow field is extremely difficult and that it is most helpful, 
if not necessary, to have a physical model in mind in order to interpret those 
measurements. In planning the experiments, as well as interpreting the results, 
the existence of large scale structures in the approximate form of interacting 
vortex rings was postulated and a consistency or contradiction with this picture 
was continually sought. 

In attempting to characterize these structures, one would first of all like 
to obtain information about their average convection velocity and their linear 
dimension. Of the two, the first quantity is by far more easily obtained with- 
out ambiguity. 

Average Convection Velocity .- In order to minimize the effect of the fine 
scale "eddies" always present within the shear layer, measurements were made in 
the potential flow surrounding the turbulent layer, where the "footprint" of 
the large structures are expected to be most prevalent. The pressure cross- 
correlation measurements shown in Fig. 3 provide the required information. The 
probe separation distance divided by the passage time, ^ , of the turbulent 
structures moving from one probe to the other is the convection velocity. Here, 
^ , is identified with the delay time corresponding to the peak correlation 
amplitude in the figure. The distribution of the convection velocity normalized 
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by the jet exit velocity along the jet is shown in Fig. 12. Along the potential 
core the ratio is in good agreement with previous measurements (refs. 3, 16 and 
17). Downstream of the potential core the convection velocity remains approximately 
one-half of the centerline velocity and therefore, decays asymptotically as 1/x . 

Length Scales Associated with the Structures .- Much more difficulty was 
encountered in attempting to obtain information concerning the average separation 
distance between neighboring structures and about their size. This is primarily 
due to the fact that according to the visual observations these structures co- 
alesce randomly in space as they move downstream; this prevents one from estab- 
lishing a reliable phase reference that is necessary to separate spatial and 
temporal variation of the structures exhibited by the velocity or pressure 
signatures. Nevertheless, some progress has been made in this direction. 

Referring to the energy spectra of the velocity fluctuations obtained at the 
two outer edges of a shear layer (ref. 4), the energy peaks are found to corre- 
spond to the passage frequency of two a 1 ready- combined vortices while the spectrum 
near the layer center exhibits a second peak at twice this frequency attributed 
to the passage frequency of the individual vortices. With this information in 
mind one would expect that in the near potential field of the jet where velocity 
and pressure fluctuations are closely correlated (refs. 3 and 17) the pressure 
spectra or the equivalent auto-correlations would give information about the time 
of passage between successive structures. The time delay to the first minimum 
of the auto-correlations of Fig. 2 is established as a measure of the average 
passage time between successive structures at various stations downstream of the 
nozzle (Fig. 13). With this information and from the knowledge of the convection 
velocity, an average separation distance, X, , between the structures can be 
calculated (Fig. 14). Separation X increases linearly with downstream distance. 

One should note here parenthetically, that it would be erroneous to identify 
the length X with the size of the "acoustic" sources. Referring to the visual 
observations and to Powell's vortex source model (ref. 7)> it is more reasonable 
to think of regions of high correlation (and high vorticity) as the coherent 
volume acting as the "acoustic" source. The dimension of this region is unlikely 
to be larger than the local shear layer thickness, and therefore, less than A . 
Unfortunately, no successful technique was found with pressure or velocity 
transducers for a quantitative measurement of this length scale; perhaps, the 
use of a vorticity meter might prove to be more profitable. 

Similarly, no additional information was provided by the pressure cross- 
correlations measured radially across the jet at positions lying on a cone angle 
of 10® (Fig. 4). The correlation is decreasing rapidly with axial distance 
from the nozzle and is essentially zero at about x/D =5. in qualitative agree- 
ment with the results of Mol lo-Chr istensen (ref. 18). In retrospect this is not 
surprising, since two pairing vortex rings are unlikely to retain axisymmetry; 
higher order modes are expected to develop in the turbulent region. 

Coalescence and Shear Layer Growth .- Interaction between neighboring vor- 
tices referred to as coalescence has been postulated to be the major mechanism 
of shear layer growth (ref. 4). Subsequently, Laufer, Kaplan and Chu (ref. 19) 
conjectured that this Interaction may be an essential contributor to the noise 
producing process. It is, therefore, of considerable importance to study this 
phenomenon . 
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The occurrence of the first pairing is demonstrated in the set of auto- 
correlations shown in Fig. 15 i measured by placing a hot wire into the shear 
layer at the radial positions where the velocity fluctuations are maximum, a 
position which should select the peak of the eigenmode. At axial stations 
0<x/D<.5 the auto-correlation looks like a damped oscillation with a Strouhal 
number of 0.23 based on the initial shear layer thickness of O.O 5 diameters 
and on half of the jet velocity. This frequency corresponds to the eigenf requency 
of the preferentially amplified, two dimensional shear layer instability. The 
implication is that the jet shear layer is essentially two dimensional near the 
exit plane, and the subsequent two dimensional instability grows with the shear 
layer until curvature effects force the local instabilities to phase lock into 
a growing axisymmetric vortex sheet. Further downstream, at x/D >.7 (Fig. 15)> 
every other wave has disappeared. This is the completion of the first pairing 
and occurs after the transition of the shear layer to turbulent flow. It should 
be emphasized that this production of an exact subharmonic is qualitatively 
different from what would be observed if dissipation or diffusion were the 
important physical mechanism here. Both of these processes lack the necessary 
phase reference to produce the discontinuous spectral shift observed in Fig. 15- 
The continuous length scale spreading produced by turbulent diffusion would 
result in a continuous spectral shift, which Is, in fact, observed beyond one 
diameter downstream. Consequently, beyond the first pairing, auto-correlation 
measurements are insufficient to distinguish between diffusion and pairing as 
the important mechanism in shear layer spreading. 

In order to determine the phase relationship across the mean coherent 
structure at positions further downstream, two-point education, described in the 
section on Special Statistical Techniques was used. The educted waveforms are 
shown in Fig. 5. Close to the nozzle (Fig. 5a) there is an abrupt 180° phase 
shift in the center of the shear layer where the fluctuation amplitudes are small. 
This behaviour is suggestive of a simple vortex sheet. Another 0.1 diameter 
downstream, however, this behaviour has changed; the phase shift is spread out 
over 0.05 cm, and in the center there are twice as many crests (Fig. 5b). 

This result is consistent with the findings in the shear layer (ref. 4) where 
the double frequency is related to the pairing process. The double structure 
is barely discernable at x/D ” 0.9 (Fig. 5c); thereafter, the phase shift 
appears smooth and continuous (Figs. 5d and 5e). 

The failure of both of the above techniques to detect the coalescence at 
downstream positions larger than a diameter is believed to be primarily due to 
the phase jitter in pairing location. In order to get an estimate of the mag- 
nitude or standard deviation of the fluctuations in pairing location the following 
calculation has been carried out based on a simple model. 

If we consider the coherent shear layer structures to be vortex rings, then 
the rate of vortex pairings within a given interval (x-j , X 2 ) is the rate n 

at which waves disappear in that interval. 
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Here f is the vortex passage frequency (based on auto-correlation period, 
Fig. 13 ). From that figure, frequency decays as 


along the first four diameters and 



If ^2 - is the mean distance, 1 (x) , for one pairing to take place, then 
half the waves passing station will have disappeared before reaching 

station X 2 : 


so that 


_/ 



- 





(9) 


and 1 (x) = X . Therefore If Is the mean location of the first pairing 

then the location of the n^^ pairing will be 








a 




( 10 ) 


A similar analysis predicts mean pairing locations beyond the potential core. 
If the first pairing were located at x^ = 0.24 diameters, then the deduced 
pairing locations would be as follows; 

n 1 2 3 4 5 § Z_ 

x/0 .24 .DO 1.3 4.4 6.3 8.9 12.6 


Downstream of the jet core the pairings become Increasingly sparse and^ of the 
seven deduced pairing locations, only three (numbers 4 through 6) occur in the 
principal noise-producing regions of the jet (ref, 8). Then, according to this 
mode), the majority oF the noise would be produced by the fourth, fifth and 
s i xth pa i rings . 

Vortex pairing has been postulated to be the major mechanism of shear layer 
growth (ref. 4), The fact that jet shear layers are observed to grow linearly 
with X should permit an estimate of the upper limit to the phase jitter assoc- 
iated with pairing* If the location of the pairing is assumed to follow 

a Gaussian distribution about its mean position, x , the shear layer thickness, 
h(x) , can be modeled as ^ 


12 



( 11 ) 



= Uncertainty in the pairing location expressed 
as a fraction of the local vortex spacing. 

Equation (11) has been computed for the jet core region for three values of 
<^7" (fig. I6a, b and c) . The shear layer spread does not become smooth until 
is greater than 0.8; however, for 1.5 the growth rate is slower 

than linear. On the basis of this model one would estimate the dispersion in 
pairing location to be of the order of the vortex spacing. This relatively 
large Jitter can indeed explain the reason for the difficulties in studying the 
pairings by conventional statistical techniques. 


Preliminary Results Using Markov Transitions 


The most controversial feature of the two-point eductions is the selection 
of a relevant trigger condition. The simple condition of level and slope is 
sufficient for detecting mean structure, whereas detection of an intermittent 
event without a fixed phase relationship to the mean structure requires subtlety. 
A promising method, derived from pattern recognition techniques, is to study the 
appropriate Joint probability distribution and transition operator in hopes of 
identifying some self-similar, deterministic subprocess associated with the 
turbulent structure of interest. Identification of the subprocess provides the 
necessary trigger condition for educing the structure with one or more signal 
probes separated in space from the trigger location. We postulate that for high 
Reynolds numbers the state of the fluid is completely described by the field 
u()<, t) , and that at a fixed position the future state of the fluid is described 
statistically by the present state; i.e. the turbulence may be represented as a 
model Markov process with a vector Langev in Equation 

The vector jf is the time-averaged transition - u{t + rTt) - u(t) given 

the state (u^ , u^). It is referred to as the expected value of the transition 
and represented formally by 

(13) 
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The vector 
That is 



is purely random with zero mean and unity variance. 






and 



(14b) 


The tensor ^ is a matrix of influence coefficients and when contracted 
against i/A represents the Brownian motion of the fluid at a location u , 
and aij^'s the familiar Kroenecker delta. Taken together, 
represents the diffusion of the system state around a location u during a 
time step i/t , and _f efx. represents the corresponding drift to the new 
average state. 

It is useful here to present some statistics of the velocity components 
u and u as discussed in the part on Joint Probability Measurements. 

Tfiei r joint probability distribution taken at two points across the shear layer 
are given in Fig. 6. The distribution at r/D = .5& (Fig. 6a) is skewed 
towards positive velocity fluctuations u and u , indicating entrapment of 

high momentum, irrotational fluid from the core, fn the third quadrant (neg- 
ative u and u fluctuations) there is a well defined "edge" to the dis- 
tribution which appears topographically like a "forbidden zone": There are 

a negligibly small number of occurrences of velocities with (u , u^) values 
lower than this edge value, and the isoprobability contours are closely spaced 
near this edge. Hence the fluid velocity spends a good deal of time next to 
the edge but is rarely able to accelerate beyond it. The distribution measured 
closer to the potential core (r/D = .49, Fig. 6b) does not exhibit such an 
edge, but shows a general consistency in that it Is skewed towards the negative 
u , u„ fluctuations. This distribution is dominated by entrainment of low 

X ’ ♦ * 

momentum fluid from the outside of the jet. 


The correlation coefficient R is tabulated in Table 1, along with r.m.s. 
values of and fluctuations. In Table 1, our results are compared 

against the tow speed 'jet measurements of ref. 3, the jet measurements of 
Bradshaw, Ferriss and Johnson (ref. 20), and the two-dimensional shear layer 
measurements of Jones (ref. 21) and Liepmann and Laufer (ref. 22). There is a 
surprisingly broad spread among published values, but generally speaking our 
measured r.m.s. levels of fluctuating velocity components is in acceptable 
agreement while our correlation coefficients are slightly low. 


The actual measurement of Joint probability densities and state space 
transitions is a difficult process with analog equipment; however, these are 
among the easiest calculations to program for digital computers. In practice, 
records of u and u are sampled periodically and stored in parallel. 

Each occurrence of ^ in a box (u ^ Au 1 1 \ u t Au /Z ) generates a 
unique index. By counting the number of such occurrences, the joint probability 
distribution is generated. At the same time one can look forward one step in 
time to produce the instantaneous differential {Au , ^u ), which is 


14 



associated with the state-space location (u , u ) . The ensemble average 
of , determined at each state-space location jj , is just the transition 

operator fdt . Higher order transition operators are measured in similar manner. 

Two sets of measured transition vectors obtained at the same two Indicated 

points in the shear layer, are shown in Figs. 17a and b. The transitions are 
indicated superimposed on isoprobab 1 1 tty contours. In Fig. 17b generally they 
point radially towards the probability peak; i.e., the most probable transition 
will be towards the most probable state. However, in Fig. 17a, near the edge 
of the probability distribution, the transition vectors (when laid end-to-end) 
indicate a swirl -like pattern. It would be useful to develop a quantitative 
measure of such subprocesses. There are two that come readily to mind; the 
inner product u • _f and the cross product ux_f . The latter would be par- 
ticularly interesting since the integral 

wouid provide a distributive measure of how much "swirl" is present. However, 
the quantity i£xf has so far eluded our physical interpretation. The Inter- 
pretation of u*£ is much more straightforward, and is a direct consequence 
of the development of the study of joint transitions. 


The Brownian motion of the system provides the leading term in the expec- 
tation value of the joint transition tensor Ju. 

^/A-fuj)u] -- 

Note that Equations 13 and 15 are the LIndeborg conditions of a Markov process. 
It is a matter of simple algebra to show that the expected transition in the 
cross product u. u. is 

where the tensor ^ in the matrix product 


a] ■ S,y - Ktj “j-f. 


06 ) 




07 ) 


J 

If Equation 16 is rewritten in terms of products of fluctuating components 
u. , 'u. and contracted over the indices, then 

is the expected transition in turbulent kinetic energy. Note that u, f is 
the inner product discussed earlier. The contracted tensor ^ ^ 


( 18 ) 




( 19 ) 


can be measured directly, and Is insensitive to the coordinate system mis- 
al ignment. 
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It is interesting to compare turbulent kinetic energy (T.K.Ej transitions 
with T.K.E. density probability distribution Itself defined as (u^ + p(u)^ 

Fig. 18 shows isoprobability contours of T.K.E, density at r/0 = .56, while 
Fig. 19 displays contours of equal values of T.K.E, transition (16) at the same 
radial position. Positive values of /“(T.K.E.) indicate transitions to levels 
of higher T.K.E. while the cross-hatched region denotes negative transitions. 

The distribution of T.K.E. transition levels is consistent with intuition in 
that regions of small T.K.E, (i.e. near the mean) show the largest positive 
Jr(T.K,E.) and vice versa. In comparing Figs. I 8 and 19> it is apparent that 
the T.K.E. peak associated with entrapment of high momentum fluid (positive 
u and u ) occurs at a state corresponding to small rates of change of T.K.E. 

By implication, this fluid spends most of its time traveling between states of 
large positive and negative acceleration, which is what one would expect. 

However, the T.K.E. peak associated with low momentum fluctuations which exhibits 
a sharp "edge” like the P(u) distribution, is located at a state where T.K.E. 
transition levels are large in magnitude. Apparently when a particle reaches 
the edge of the energy distribution, it experiences a large acceleration towards 
lower energy levels. This acceleration must be parallel to the edge and is 
consistent with the shape of the P(u) ridge shown in Fig. 6a as a dotted line. 

The set of data P(u) , X(h) .^^(u) provides a detailed and sensitive 

indicator to the state of the turbulent fluid. in contrast, both time correla- 
tions and spectra taken at the two locations cited show little difference or 
distinctive features. Both of these quantities can be extracted from F, 
and ^ , although with some difficulty. 

The advantage to these measurements is that they provide inputs to newer 
theories of turbulent structure involving Markov processes, and hence can prove 
to be useful to the analyst. 


Artificial Initial Disturbances 


The experiments of Crow and Champagne (ref. 2) have shown that the gross 
properties of the jet, such as axial velocity distribution and spreading rate, 
can be influenced by controlling the initial conditions of the jet. Their 
forcing technique— acoustically driven stagnation chamber resonance - requires 
relatively high disturbance amplitudes, typically in excess of 1% of the jet 
exit velocity. The method adopted in the present work consisted of generating 
an annular shaped acoustic wavefront (Fig. 7) that surrounds the entrainment 
region of the jet and that has a wave number vector parallel to the jet axis. 

It was hoped that by driving the shear layer in this way, at the right frequency, 
a given vortex ring pairing could be phase locked with the disturbance, thus, 
facilitating detection. 

Results have shown that the above described forcing technique apparently 
did inhibit the phase jitter but only during the first and possibly the second 
vortex pairing, by phase locking the laminar eigenmode. As a consequence, the 
transition to turbulence was delayed from 0.3 to 0.5 diameters. 
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Although this type of driving has no obvious influence over mean proper- 
ties of the jet (such as axial and radial velocity profiles) there is an 
interesting influence over shear layer turbulence. Broadband turbulence levels 
in the transition region are increased by forcing, probably by strengthening 
the eigenmode. Further downstream, however, forcing is observed to decrease 
considerably the shear layer turbulence. Fig. 20 shows radial traverses of 
the forced and unforced shear layer at the downstream station x/D = 4. Forc- 
ing was at 3*6 kHz, (the laminar eigenmode) and the mean velocity profile is 
shown for orientation. The turbulence level is reduced by 50% at the center- 
line and by 20% at the peak value. It is to be noted that this result is In 
direct contrast to the effect observed by Crow and Champagne (Fig. 28, ref. 2), 
This observation emphasized that driving the jet core produces qualitative 
differences from driving the shear layer turbulence only, 

It-was, of course, of considerable interest to investigate whether or not 
the reduction in turbulence level effects the far pressure field. The above 
method of excitation was repeated, therefore, in a jet placed in the anechoic 
chamber. At a Mach number of 0.5 intensity and spectral measurements indicated 
no discernable change when the artificial disturbance was introduced into the 
shear layer. Instead of continuing the far field investigation, it was decided 
to attempt to understand in more detail the role played by the artificial dis- 
turbances in the near field. 


tnter-Relat i onsh i p Between the Large Structures and the Noise Field 


As pointed out in the Introduction, the question of the large scale 
structures as potential acoustic sources is a relevant one. A direct experi- 
mental confirmation of this problem is a most difficult one; the evolution in 
time of large coherent flow regions would have to be measured, (simultaneously 
with the far field pressure) a task that present techniques cannot as yet 
perform. A more indirect approach to the problem has been adopted in the 
present work: one postulates a feasible noise mechanism associated with Lhe 

large structures and one makes certain measurements in the far field to look 
for the consequences of the assumed mechanism. This requires a highly direc- 
tional microphone In order to be able to detect the noise generated over a 
relatively small region of the jet. 

Assuming that one of the dominant mechanisms of the jet noise Is associated 
with the coalescing vortical structures in the flow, one would expect to find 
certain features in the far field. In particular, one can think of at least 
three such character! st ic features: 

1. the noise produced by any single region of the jet ought to be 
intermittent; 

2. evidence of spatial coherence should be present in the noise 
field, since the coalescence is a spatially coherent process; 

3. the noise spectrum should reflect some characteristic time that 
is associated with the coalescence process. 
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In studying the data obtained with the directional microphone with re- 
spect to these features the following comments can be made. At subsonic jet 
velocities no i ntermi ttency has been observed in the far field; it is quite 
possible that at these velocities the intermittency is being obscured by the 
spatial averaging of the mirror's finite focusing window. It should be men- 
tioned, however, that at supersonic velocities highly intermittent signals have 
been observed which in fact are believed to emanate from the large structures. 

With respect to point 2, spatial coherence measurements have been made in 
the far field with a nine microphone array located at equally spaced intervals 
along the 36 in. line indicated in Fig. 21. Spatial coherence is expressed as 
cross-correlation coefficients with respect to the center microphone. Polar 
angles of 30° and 90° were chosen. Although moderate coherence has been 
found the evidence was judged to be inconclusive. 

A more interesting result has been found related to the characteristic 
time of the coalescence. According to the model of Winant and Browand (ref. h) , 
the pairing time is proportional to the vortex size (which in turn is propor- 
tional to the vortex spacing) and inversely proportional to the convection speed. 
But the ratio of vortex separation distance to convection speed corresponds to 
the average passage time between successive vortices (Fig. 13). These data 
have been replotted in the form of a non-dimensional passage frequency in Fig. 

22 and compared to the following far field measurements. Narrow band noise 
intensity distribution measurements were made with the reflector microphone 
located 90° to the jet axis and directed toward narrow segments of the jet 
along its axis. Fig. 23 shows a set of intensity distributions as a function 
of x/D . It is seen that at a given x station one particular frequency (or 
Strouhal number) contains more energy than any other. These non-dimensional 
frequencies have been plotted as a function of x/D in Fig. 22. It is seen 
that the near and far field data show surprisingly similar trends. (Unfortu- 
nately, the accuracy of the far field measurements near the jet exit is limited 
by the spatial resolution problems.) This implies that the time scales associ- 
ated with the noise radiated normally from a given region of the jet is related 
to a time scale associated with the near field temporal coherence. Clearly, 
such a conjecture is insufficient to provide a convincing argument to the above 
model and much more work is necessary to clarify this question. 

CONCLUSIONS 


Various statistical measurements carried out within and near a turbulent 
jet are shown to be consistent with the conjecture that large scale organized 
structures do exist in such a shear layer; they interact with each other con- 
sistent with visual observations made at low Reynolds numbers. Unfortunately, 
because of their random interaction by coalescence, only a limited amount of 
information about their nature and behaviour could be extracted even with the 
most modern statistical techniques. 

In particular, the average separation distance between neighboring 
structures was obtained from conventional auto-and cross-correlation measure- 
ments. It was found to increase linearly with x along the jet axis. 
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The pairing process, clearly observable by visual techniques, could be 
detected only in its earliest stages; more sophisticated statistical methods, 
such as two-point eduction or state-space statistics were unable to identify 
the occurrence of coalescence further downstream in the flow. It is argued 
that the difficulty is primarily due to the random occurrence of the coalescence 
which prevents one from acquiring a reliable phase reference necessary for the 
detection. 

The presence of the structures and the important consequences of their 
interaction was further brought out by artificial excitation of the jet. Al- 
though no measurable changes were detected in the mean velocity field of the 
jet, a considerable decrease in the velocity fluctuations could be attained at 
four diameters downstream of the jet. It is conjectured that a particular 
eigenmode of the shear layer can inhibit the interaction process between the 
vortex rings and thereby delay transition. 

Finally, it should be noted that as yet the work has not yielded a tech- 
nique that would allow one to study directly the connection, if any, between 
the large structure behaviour and the far field noise. The only circumstantial 
evidence produced by the present research is the observation that the passage 
frequency {which is believed to be proportional to the rate of subharmonic 
formation) of the structures varies with x in a surprisingly similar manner 
as the frequency corresponding to the maximum intensity radiation emanating 
from the same value of x. It is evident that this point, a most important one 
for understanding the jet noise problem, needs to be investigated in much more 
detail. 
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h(x) 

k 

M(x) 
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M 
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p(x) 

p' 

R 
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U,(x, t) 
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9 
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U(x) 

U 

e 

X 
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son?c velocity 

sonic velocity - ambient 

jet exit diameter 

Jet exit radius 

wave frequency 

transition vector (expectation of velocity partial time derivative) 
Jet shear layer thickness 
wave number 
local Mach number 

local Mach number - referred to ambient 

exit plane Mach number 

joint probability distribution 

static pressure fluctuation 

r.m.s* pressure fluctuation 

correlation u u / u* u* 
x r X r 

Reynolds number U 0/t-^ 

€ 

near field time scale 
local velocity (tensor notation u, ) 
local velocity axial component 
local velocity ^ radial component 
fluctuating local velocity ^ 

X 

r.m.s, velocity fluctuation = (u.^-* u.^)^ 

average turbulence convection velocity 

phase velocity 

jet centerline velocity 

jet exit plane velocity 

field position vector 



X 


- axial position 


/q 

e 

X (x) 



'T 

0 


= radial position 

= difference in radial position from shear layer centerline 

= d{t + ^t) “ Q.{t) for any Q. 

state space joint transition tensor 

spherical polar angle 

local wave length 

kinematic viscosity 

local density 

dimension uncertainty in vortex pairing location (Eq. 11) 
expectation of square of velocity time derivative (tensor quantity 

•J 

time delay 
azimuthal angle 
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TABLE 1 

Comparison of Turbulence Levels 


Jet Measurements 
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27 
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Bradshaw, Ferr i ss , Johnson 
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Lau, Fisher, Fuchs 
(Ref. 3) 
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2-Dlmens Iona) Shear 
Layer Measurements 








Jones (Ref. 21) 
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. 18 
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Liepmann, Laufer 
(Ref. 22) 

18.5 

.26 

. 16 

.54 

.17 

.12 

.56 




Auto^corre lat ion , r — , (dimensionless) 



Time delay, , msec 
(a) Axial stations, x/D, 1,0 to 4.0. 



Figure 2. -Auto-correlation of pressure signals measured at various stations 
along edge of jet. The arrows (^ ) denote the time of the first minimum. 
Radial positions lie along 10° cone. 
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(b) Axial stations, x/D , 6 to 10 
Figure 2. -Concluded. 





Cross Correlation, PaPd(TTK (arbitrary units) 



Time delay, “ZT, msec 
(a) Axial Stations, x/D , .5 to 1.5, 



Figure 3* "Axial cross-correlations of pressure signals at various stations 
along edge of jet. Arrows (f) denote time of maximum cross correlation. 
Radial positions lie along 10° cone. 
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Figure 3* " Concluded, 
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Figure 4. 
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- Maximum of azimuthal cross-correlations of pressure 
measured across the jet along a cone angle of 10°. 
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(b) x/D = .28 ; ^ r/h , -O.I3 to 0.07 

Figure 5. - Continued* 




Educted Velocity, u , (arbitrary units) 
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(c) x/D = 0.9 ; 4 r/h, -0.16 to 0.63 

Figure 5. - Continued. 




Educted Velocity, u , (arbitrary units) 



Relative Time (msec) 

(d) x/D = 4 ; r/h, -0.23 to 0.46 
Figure 5. ^ Continued* 




Educted Velocity, u , (arbitrary units) 








Radial velocity fluctuation, u^ - , Radial velocity fluctuation 

(relative units) (relative units) 



Axial velocity fluctuation, u^ > u^, (relative units) 
(a) r/D = 0.56; u/U = 0.60 



Figure 6. - Joint Probability Distributions P(u , u ). Isoprobability Contours 
shown for x/D = 2. Contour level shown as fraction of the distribution peak. 










Radial position, r/D, Centerline velocity, 

(dimens ion less) (dimens ion less) 






Axial Velocity, u /U , (dimensionless) 



Radial position, r/D , (dimensionless) 
Figure 10, - Radial Velocity Profiles. 
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Acoustic frequency, (dimensionless) 

Figure 11, - Dimensionless overall power spectral density of jet noise* 
Mg 0*90 + 0*02, Reynolds numbers 0,88 to 3,3 x 10^. 
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Axial position, x/D , (dimensionless) 
Figure 12. - Turbulence Convection Velocity Decay. 
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Figure 13. - Auto-Correlation Time Scale (T) vs. Downstream Station 
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Downstream station, x/D , (dimensionless) 
Figure 14. - Spacing Between Coherent Structures, 
a Ko & Davis (Ref. 16) 

0 Fuchs (Ref. 17) 

V Lau, Fisher, 6- Fuchs (Ref. 3) 
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Figure 16. - Concluded. 




Radial velocity fluctuations, u - u , (relative units) 
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Figure 17. - Experimentally measured expectation of velocity derivative f 

(Eq. 13). The relative magnitude and direction of jF are indicated 
by the arrows (f), and the pattern is modulated by” the isoprobability 
contours (fig. 6). 



Axial velocity fluctuations, u - u 
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Axial velocity fluctuations, u-Z , (relative units) 
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(fa) r/D = 0.49; x/D =2.0 


Figure 17. - Concluded 





Axial velocity fluctuations, u - u , (relative units) 
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r/D = .56; x/D = 2.0 

Figure 18. - Turbulent Energy Isoprobability Density Contours. 



Figure 19. - Turbulent Energy Transition Levels. Cross-Hatched Region Indicates 
Negative Transition Levels. Note: Contours at Equal Intervals. 
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Radial position, r/D, (dimensionless) 

Figure 20. - Effect of Forcing on Shear Layer Turbulent Profiles, x/0 
Mean Profile included for Reference 
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R.M.S. fluctuating velocity, u*/U , (dimensionless) 




Polar Angle, degrees 

Figure 2t. - Far Field Acoustic Spatial Coherence. 
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Figure 23. - 
along a I" 
far field 
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Axial location, x/D , (dimensionless) 

Processed, narrow band acoustic source intensity distributions 

» = .97 jot, as seen by a reflector microphone in the 

normal to the jet axis. 



